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ABSTRACT 

We assess the possibility of determining the Hubble constant Hq by measur- 
ing time delays between multiple images of supernovae gravitationally lensed by 
rich clusters of galaxies and combining these delay measurements with detailed 
cluster-potential models based on other lensing constraints. Such a lensing 
determination of Hq would be complementary to those obtained from galaxy- 
QSO lensing studies, and could potentially be better calibrated. We show that 
relatively low-redshift {z ~ 0.2), significantly elliptical clusters have appreciable 
lensing cross sections for observable image pairings with tractable time delays on 
the order of a few years despite large lensing mass scales. We find that a tar- 
geted search for such image pairs would be a significant undertaking for current 
observatories, but that it would be appropriate for a facility such as the proposed 
Large-aperture Synoptic Survey Telescope. 

Subject headings: distance scale — gravitational lensing — galaxies: clusters: general- 
sup ernovae: general 



1. Introduction 

The best measurements of the Hubble constant {Hq) based on the Cepheid-calibrated 
distance ladder and on the cosmic microwave background agree impressively with one an- 
other: Hq = 72 ± 8kms~^Mpc~^ and Hq = 72 ± 5kms~^Mpc~^, respectively (Freedman 
et al. 2001; Spergel et al. 2003). However, these values are inconsistent with much lower 
values of Hq based upon measurements of gravitational lens time delays and the assumption 
that lensing galaxies have extended dark halos (Kochanek 2002a,b, 2003). In this paper we 
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explore the possibility of determining Ho by monitoring known strong-lensing galaxy clusters 
and measuring time delays between multiple images of high-redshift supernovae (SNe). 

The number of known giant-arc clusters is encouraging. Wu et al. (1998) present a 
summary of 38 strongly lensing clusters from the literature, containing a total of 48 giant 
arcs and arclets (see also WiUiams, Navarro, & Bartelmann (1999)). More recently, Luppino 

et al. (1999) have reported the discovery of 8 more giant-arc clusters, Gladders, Yee, & 
Ellingson (2002) 6 more, and Zaritsky & Gonzalez (2003) 3 more. Recent observations of the 
galaxy cluster Abell 1689 with the Advanced Camera for Surveys (ACS) on the Hubble Space 
Telescope have revealed at least 30 new multiply-imaged background sources, with an average 
of ~3 images each (D. Coe, 9 January 2003 presentation "Deep ACS and Keck Observations 
of A1689", in conference "Gravitational Lensing: a Unique Tool for Cosmology", Aussois, 
Savoie, Prance). All these image systems (together with previously known strongly lensed 
features) provide constraints on the projected gravitational potential of the cluster — many 
more constraints than a quadruple-image quasar lens system places on its lensing galaxy. In 
addition, since multiple strongly lensed sources behind a single galaxy cluster in general lie at 
different redshifts, the cluster potential map will not suffer the "mass sheet degeneracy" that 
afflicts galaxy-quasar lens systems (Saha 2000). Although A1689 is perhaps the strongest 
cluster lens in the sky, we may assume that other known strong-lensing clusters will likewise 
exhibit many new multiple-image systems when observed with ACS and later instruments, 
and detailed cluster modeling like that of AbdelSalam, Saha, & Williams (1998a,b) and 
Kneib et al. (1995, 1996) but with increased accuracy will become possible. Thus although 
the structure of galaxy clusters is in general more complex than that of individual galaxies, 
the projected potentials of clusters may in the near future be more tightly constrained by 
lensing than those of individual galaxies, and cluster-lens time delay would therefore be more 
precisely calibrated. (See Schechter (2000) for a general discussion of the determination of 
Hq from gravitational- lens time delays.) 

This work is distinguished from previous studies by its quantitative focus on the pos- 
sibility of measuring SN-image time delays in cluster-scale lens systems, and combining 
them with additional strong-lensing constraints to determine Hq. The original proposition 
to measure Hq using gravitational lens time delays envisioned distant supernovae (SNe) as 
the variable sources to be lensed by intervening galaxies (Refsdal 1964); quasars and their 
variability had yet to be discovered. Kovner & Paczyhski (1988) examined the possibility 
of observing SNe exploding within giant-arc source galaxies. They made a detailed analysis 
for a SN very near to and inside a caustic "cusp" and showed that for such configurations 
the time delay between multiple images could be as short as days to weeks, although they 
did not explicitly consider the possibility of making an Hq determination. Giraud (1992) 
reported the detection of variability in two separate arclets of the cluster CI 0302-1-1658, 
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consistent with a single source and a time delay of 1 to 2 years. Various papers have exam- 
ined magnification and time delay effects for SNe in field surveys lensed by intervening halos 
across a spectrum of mass scales (Linder, Schneider, & Wagoner 1988; Marri & Ferrara 1998; 
Dahlen & Fransson 1999; Marri, Ferrara, & Pozetti 2000; Porciani & Madau 2000; Holz 2001; 
Goobar et al. 2002; Oguri, Suto, & Turner 2003; Oguri & Kawano 2003), and several have 
explicitly considered the possibility of constraining Hq using SN time delays. However, all 
these works have modeled the lenses as circular, and in this paper we show that circular mod- 
els seriously under-predict the capacity for cluster-scale lenses to produce multiple images 
with reasonably short time delays (as one might anticipate based on the results of Kovner 
& Paczyhski (1988) for sources very near to the caustic cusp of a lens). Other studies have 
considered targeted searches for lensed SNe in galaxy cluster fields (Miralda-Escude & Rees 
1997; Kolatt & Bartelmann 1998; Saini, Raychaudhury, & Shchckinov 2000; Sullivan et al. 
2000; Gal- Yam, Maoz, & Sharon 2002; Gunnarsson & Goobar 2002), but their quantitative 
focus has been on pure SN detection rates and not on the prospects for measuring time 
delays that we quantify here. With the exception of Saini et al. (2000), these works have 
also approximated clusters with circularly symmetric lens models. 

This paper is organized as follows. Section 2 presents a brief discussion of relevant 
gravitational lens theory and defines the particular lensing cross section that is central to 
our calculations. In Section 3, we examine the lensing behavior of several cluster models 
in the context of observational constraints on minimum magnification and maximum time 
delay. Section 4 discusses the accuracy in the determination of Hq that we might hope 
to obtain from individual cluster time delays, both for a simplified model and for realistic 
models. Section 5 presents a calculation of detection rates for observable SN image pairings 
with acceptable time delays in the strong-lensing region of a specific cluster. Our approach in 
this section is similar to that of Sullivan et al. (2000) and Gunnarsson & Goobar (2002) (and 
differs from that of Saini et al. (2000)) in that wc do not directly consider the possibility of 
SNe within known lensed galaxies, but rather attempt a prediction of SN detection rates in a 
targeted cluster field based on based on an assumed cluster mass model and cosmic supernova 
rate density function. (Lanzetta et al. (2002) argue that above redshifts of ~1.5-2, even the 
Hubble Deep Field observations are insensitive to the surface brightness of most rest-frame 
ultraviolet emission.) Finally, in Section 6 we discuss the observational outlook in view of 
the calculations of Section 5. We find that a measurement of Hq based on cluster time delays 
would be challenging but not impossible with today's telescopes and instruments, and that it 
would be an appropriate project for a large telescope operating in a dedicated survey mode 
such as the proposed Large-aperture Synoptic Survey Telescope (LSST) (Tyson et al. 2002). 

Throughout this paper, we assume a fiat, vacuum-dominated cold dark matter (ACDM) 
cosmology with matter and vacuum density parameters Qu — 0.3, fl\ — 0.7. For the Hubble 
constant, we take Hq = 70 /lyokms"^ Mpc~^ with h-jo = 1. 
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2. Lensing Theory 

Gravitational lensing can be described elegantly through the application of Fermat's 
principle (Schneider 1985; Blandford & Narayan 1986). For a source at (2- vector) angular 
position P, the time of arrival relative to the unlensed case for image positions 9 is given by 



^ ^ (1 + ^l) DlDs 
c -Dls 



(1) 



where zj^ is the lens redshift, D]^, D^, and £)ls ^-re angular diameter distances to the lens, 

— * 

to the source, and from the lens to the source, and tjj{6) is proportional to the Newtonian 
potential of the lens projected perpendicular to the line of sight (Narayan & Bartelmann 
1996). Images of the source will be seen at positions 6 where the arrival time is stationary: 
that is, at solutions to the "lens equation" obtained by setting the gradient with respect to 
9 of (1) to zero: 

^{e) = e-VfP{0) . (2) 

The scalar magnification /i (the flux ratio of lensed to unlensed images of a source) is given 
by the ratio of lensed to unlensed differential angular area, thus it is given by the Jacobian 
determinant of the mapping 6(/3) (the local inverse of (2)): 



fj, — det ( 



fdOi 



det 



(3) 



For an excellent presentation of gravitational lensing in more detail, we recommend the work 
of Narayan & Bartelmann (1996). 

The detection of multiply imaged SNe will be limited by both attainable photometric 
depth and maximum tolerable time delay. Accordingly, we define a{iija.m, ^^max) to be the 
angular cross section in the source plane for lensing into a pair of images with the fainter of 
the two having (absolute- value) magnification of at least /Xmin and the time delay between 
the two being at most Atmax- For a particular physical lens model, this cross section will 
depend upon the redshifts of the lens and the source, but we will suppress this dependence 
in our notation. 



3. Lens Models 

In this section we examine cr(/Xmin, Atmax) for three progressively more detailed lens 
models, and discuss the issue of normalization for the purpose of comparing between them. 
These models have been analyzed in detail by Kassiola & Kovner (1993) and by Kormann, 
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Schneider, and Bartelmann (1994). Particularly convenient formulas are given by Keeton 
& Kochanek (1998). We refer the reader to these papers for analytic lensing potential ex- 
pressions and for detailed discussions of the models' lensing behavior. Wherever numerical 
solution of the lens equation (2) is needed, we use our own implementation of the grid-based 
numerical lens-equation solution algorithm with 2D Newton-method solution refinement de- 
scribed by Schneider, Ehlers, & Falco (1992). 



3.1. Singuleir Isothermal Sphere 

The simplest reasonable model for an extended astrophysical mass distribution is the 
familiar singular isothermal sphere (SIS). It is parametrized solely by its velocity dispersion 
Uy. The scaled, projected mass distribution (convergence) of the SIS is given by 

Here, E is the physical surface mass density of the lens and Ecr = (c^/47rG)(Ds/-DL-DLs) 
is the so-called critical surface mass density, is the "Einstein radius" that sets the sole 
angular scale of the model for given lens and source redshifts: 



= 47r^^ (radians) = (28.8arcsec) — . (5) 

c2 Ds ^VlOOOkm/sy Ds 

The lensing potential of the SIS is simply '^{ff) = 9e\0\. 

We can obtain crifiunu- Atmax) for the SIS analytically. If the angular distance P from the 
source (i.e., SN) position to the lens center is less than 6e, two images will be observed along 
a line on the sky through the source position and the lens center: one at a distance 6e + 
from the lens center (in the direction of the source) with magnification 1 + 6-e/P and one at a 
distance 9e — P (in the direction opposite the source) with magnification 1 — P (Narayan 
& Bartelmann 1996, for example). This second, fainter image corresponds to a saddle point 
of the arrival time function (1); its negative magnification signals a reversal of image parity. 
It is this magnification that determines the //min dependence of cr(//niin, A^max) as defined. If 
we substitute the solutions ^e±/? into (1) (taking into account vectorial considerations) and 
form the difference, we find that the time delay between the two images is 

A*^?(l±iil^../3.24 . (6) 

c -Dls 9-e 

where we have defined a characteristic lensing timescale by r as implied. For a cluster velocity 
dispersion — 1000 km/s, r increases steeply at first with increasing source redshift, then 
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levels off around ~ 25 h-jQ years for 2;l = 0.1 or ~70-80 h-jQ years for zi, — 0.4. (In Section 3.2 
we will relieve our sense of discouragement over such long timescales.) 

We see that for the SIS a given /3 corresponds to a unique At as given by (6) and 
a unique (absolute- value) magnification of 9e/ P — 1 for the fainter image, so the form of 
c(/imin, A^max) is particularly simple: we have two singly limited cross sections given by 

0-(/Umin) = 7r6'|(/imin + 1)^^ (7a) 

(t(AW) = l^OliAtr^^/rf , (7b) 

and 

0'(yWmm, Atmax) = Hlin [^(/imm), O-(Atmax)] • (8) 

That is, for a given (/imin, Atmax), we are either magnification-limited or delay-limited. 
(Arguably the maximum possible double-image cross section irO^ should be enforced in (7b), 
but it is built into (7a) and thereby propagates to (8).) 



3.2. SingulEir Isothermal Ellipsoid 

By making the isodensity contours of the SIS elliptical with a minor-to-major axis ratio 
q, we get the singular isothermal ellipsoid (SIE) model: 

= ^ (q) 

As noted by Kormann et al. (1994), the total mass enclosed within a given isodensity contour 
remains constant with changing q at fixed when k is expressed in this form. We adopt 
this normalization as a basis for comparing models of differing axis ratios, with 9^ still given 
as a function of and source and lens redshifts by (5). 

The introduction of ellipticity leads to richer lensing phenomena. Most significantly, 
we acquire a cross section for quadruple imaging. For axis ratios q less than 1 but greater 
than about 0.394, four images of a source will form if the source position lies within a 
diamond-shaped "tangential" caustic surrounding the lens center and inside the original 
border between singly and doubly imaged regions (the "radial cut"). Labeling these quad 
images in order of increasing arrival time, 1 and 2 are minima of the arrival time function and 
3 and 4 are saddle points. 2 and 3 are in general of much greater absolute magnification than 
1 and 4. We now have the possibility of three independent image pairings. For q <k, 0.394, 
the major-axis cusps of the tangential caustic extend outside the radial cut and are referred 
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to as "naked cusps" ; sources within the cusps but outside the radial cut will be triply imaged. 
Kassiola & Kovner (1993) Figure l(d,e) illustrates naked and non-naked cusps. 

The SIE model has an angular scale invariance that allows us to generate the function 
c(A*min, ^^max) oucc for a givcu and rescale it as needed for any source/lens redshift 
combination. With circular symmetry broken, we would not necessarily expect the doubly 
limited cross sections for the various image pairings to be of the form (8), but our numerical 
calculations show that such a form in fact yields a very good approximation. The solid 
curves in Figures 1 and 2 show ^(/Xmin) and (T(Atmax) for the three independent quad image 
pairings of least time delay for an SIE lens model with axis ratio q = 0.65. o-(/imin, Atmax) 
can then be constructed for each image pairing as per (8). For a given q, cr(/imin, A^max) of 
any quad image pairing is limited by the full angular area of the region in the source plane 
enclosed by the tangential caustic. This area increases as the lens becomes more elliptical, 
and is approximately 0.095 9^ for q — 0.65 as seen in Figures 1 and 2. 

By comparing the sohd SIE curves of Figures 1 and 2 to the dotted SIS double-imaging 
curve, we see that that the introduction of ellipticity significantly increases the cross section 
for multiple-image lensing with relatively short time delays as compared to the circularly 
symmetric SIS-lens case, at the cost of reduced fainter- image magnification ^. Of the three 
independent image pairings in the quad case, the 2nd/3rd image pairing is most significant 
for affording both large magnification and short time delay; these are the two images that will 
merge and annihilate if the source crosses outside the tangential caustic. When the source 
approaches the center of the lens, the four images form a cross and the 2nd/3rd image time 
delay approaches its maximum value. This limit can be found by solving the lens equation 
in closed form for the special case of a source directly behind the lens center: 

Aimax, 2 to 3 = 2{q-^ - q) (^^^^^"^^ V 1 " " arctau^ \/q~^^-^ ■ (10) 

For q = 0.65, this is about 0.14r as seen in Figures 1 and 2. For a lens redshift Zi^ <~ 0.1, the 
lensing timescale r as defined in (6) will be <~ 25 hjQ years for a cluster velocity dispersion 
ay = 1000 km/s, and all 2nd/3rd image pairings of a g = 0.65 lens will have time delays of 
less than ~ 3.5 h^g years. Thus although sources very near to the caustic will indeed have 
very short time delays as shown by Kovner & Paczyhski (1988), sources anywhere within the 
tangential caustic can still be lensed into multiple images with observationally acceptable 



^Here we axe comparing the high-magnification behavior of a circular lens to that of an elliptical lens. 

For any particular elliptical lens, a quad configuration will be of greater total magnification than a double, 
since the quadruple-imaging cross section carves out the very highly magnified region of the source plane 
directly behind the lens. Hence observed quad-image QSO lenses generally have greater magnification than 
doubles. 
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time delays. Furthermore, as the 2nd/3rd image time delay approaches its maximum, the 
lst/2nd and 3rd/4th image time delays approach zero for reasons of symmetry, so (10) is in 
fact only an upper limit on the shortest time delay experienced by any quadruply imaged 
source. 

With all other parameters fixed, both r and ^| scale as the fourth power of the cluster 
velocity dispersion o;^. Therefore if the axes of Figure 2 were labeled in units of years and 
square arcseconds (instead of in units of r and 6"^) , an increase in would rescale both axes 
by the same factor. This fact together with the differing concavities of the curves shown in 
Figure 2 indicates that in delay-limited cases, increasing ay will increase an SIE-quad cross 
section, but decrease an SIS-double cross section. 



3.3. Nonsinguleir Isothermal Ellipsoid 

A singular isothermal core is perhaps an unrealistic feature to assume in a galaxy cluster. 
To investigate sub-isothermal core behavior, we can add a core of angular radius s to the 
SIE to obtain the nonsingular isothermal ellipsoid (NIE): 



2 + qei + q-Wy ■ ^ ' 

To compare the lensing behavior of models with differing core radii, we must take care to 
adopt a sensible normalization. In the limit s — > 0, 6 may be identified with O-^ from above, 
but to compare singular and nonsingular models the normalization h — O-e^s inappropriate. 
Instead we take 

h^^f^[sq+^|s^ + f^ , (12) 

where 9-^ is the Einstein radius (5) evaluated for a source at some fiducial redshift Zs- 
With this normalization, k retains the proper scaling with changing source redshift, and 
the tangential critical curve crosses the major axis of the lens at fixed position with 
varying s for a source at redshift zs. This choice is meant to respect (approximately) the 
constraint that would be placed on the cluster mass distribution by a giant arc at redshift 

Several changes occur when we increase the core radius from zero. First, the model 
now has two angular scales, h and s. For a given lensing cluster, h will change with varying 
source redshift while s remains fixed. The scale invariance of the lens is broken, and the 
lensing behavior changes qualitatively between source planes at different redshifts; we now 
need to compute cross sections for a range oi s/b and to scale them appropriately. Second, in 
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multiple-imaging configurations, an additional image that was previously infinitely demag- 
nified in the singular core becomes in principle observable (accordingly, the "radial cut" of 
the singular case is now the "radial caustic"). This new image corresponds to a maximum of 
the arrival-time function (1). In this study we exclude these "maximum" images from con- 
sideration so as to have continuity of the total cross section o"(/imin, Atmax) from the singular 
case when summing the contributions of all independent image pairings. (Actual core images 
could also get lost in the hght of a cluster cD galaxy.) Finally, although the overall cross 
section for multiple imaging decreases with increasing core radius through a contraction of 
the radial caustic, the area enclosed by the tangential caustic actually increases slightly. The 
decreased curvature of the arrival-time surface in the central regions of the lens also leads to 
larger magnifications and shorter time delays. Thus the cross section a{fj,min, A^max) at fixed 
/^min, Atinax tcuds to iucrcasc for all image pairings as the core radius is increased, until the 
formation of naked cusps begins to reduce the cross section for the 3rd/4th image pairing. 
The dashed lines in Figures 1 and 2 show the cross sections for the same image pairings as 
in the SIE case above, but for an NIE lens at the so-called "umbilic catastrophe" and for 
sources at the same redshift as the normalization redshift zs- Here the umbilic catastrophe 
refers to the particular s value at which the expanding tangential caustic completely engulfs 
the shrinking radial caustic; see Kassiola & Kovner (1993) Figure 1(c) for an illustration. 
The NIE cross sections for the lst/2nd and 2nd/3rd image pairings are increased over the 
SIE case at all values of /Xmin, A^max- All 2nd/ 3rd image pairings now have time delays 
<~ 0.06T. 



3.4. More Complicated Models 

The detailed lensing behavior of a real galaxy cluster cannot be captured entirely by 
the models considered above. The effects of unmodeled cluster substructure are likely to 
be somewhat similar to the effects of increasing cluster ellipticity relative to the circular 
case: higher image multiplicity, shorter time delays between multiple images, and reduced 
magnification. In fact, ellipticity may be thought of as the leading order of substructure 
beyond circular symmetry. By considering the perhaps exaggerated axis ratio q = 0.65 we 
hope to approximate unmodeled substructure effects. An obvious refinement would be to 
carry out our cross section calculations using a more detailed and realistic cluster mass map, 
but we defer this possibility to future investigations. Our primary goal is to estimate the 
detection rates of multiply imaged SNe in a single cluster field that would be useful for the 
determination of Hq, and thus we confine our attention to the preceding models. 
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4. Accuracy of Ho Determination 

First we will examine the usefulness of a measured time delay in determining Hq within 
the simplified context of the SIE lens model. By the convenient result published by Witt, 
Mao, & Keeton (2000) for self-similar isothermal lenses, the time delay between images A 
and B is given by 

AUB = ^^(i+.L)(i4r-i^Ar) , (13) 

with the image positions measured relative to the lens center. cHq^ factors out of the D^s on 
the right-hand side. The dependence on the cosmological parameters (I^Mi^a) is generally 
weak but worth noting. For example, taking a lens redshift of 0.1 (0.3) and a source redshift 
of 1, Di^Ds/Di^s changes only by about 1% (3%) as Oa goes from to 0.7 in fiat universes. 
(If we allow for unknown cosmological density parameters, then our time delay will not give 
us Hq, but rather the distance measure Dl-DsZ-Dls-) Assuming the image positions, source 
and lens redshifts, and time delay can be measured with high accuracy, the greatest source 

— * 

of error in an Hq determination is due to the uncertainty SOq in the position of the center of 
the lens. Prom (13), this is 

^^^(,- ^_ (14) 

For positions of 2nd/3rd image pairings obtained numerically, the dimensionless vector V 
ranges in magnitude from 2 for source positions near the lens center to almost zero for source 
positions just inside the tangential caustic. V is approximately of unit magnitude for source 
positions of area-weighted median fainter-image magnification. The fractional error in an 
Hq determination is then roughly equal to the uncertainty (parallel to the image separation 
vector) of the lens center position in units of the Einstein angular scale 9e defined in (5) 
and evaluated for the source redshift of the SN. However, the magnification of both images 
increases as the source approaches the tangential caustic, so a magnification bias favors the 
detection of events with smaller fractional uncertainties in the derived Hq. 

Of course, our real hope would be to use not an SIE model (nor even an NIE model) but 
a detailed lens model that takes full account of lensing constraints hke giant arcs and multiply 
imaged galaxies. In this more general case, the fractional uncertainty in a derived Hq value 
will be related to the fractional uncertainty in the lensing potential ip{9) and the square of 
its gradient, evaluated at the two SN image positions. (This can be seen by substituting 
(2) into (1) to eliminate the source position j3.) Based on their nonparametric lensing 
inversions, AbdelSalam et al. (1998a,b) report fractional uncertainties in the convergence 
K of < 20-25% in the strong lensing regions of interest in the clusters Abell 370 (using 6 
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multiple-image systems) and Abell 2218 (using 3 multiple-imaged systems and a number 
of singly imaged, distorted arclets). They demonstrate clearly that the mass distribution 
becomes more tightly constrained as the number of nearby lensed images increases. The 
dramatic increase in multiply imaged sources that ACS and later instruments can uncover 
should certainly push the uncertainty in reconstructed cluster mass distributions well below 
10% in the inner cluster regions. Furthermore, strong-lensing features constrain the potential 
directly, and the potential is an integral over the mass distribution. Hence we expect the 
uncertainty in the cluster potential (the more relevant uncertainty for Hq determination) 
to be less than the uncertainty in the reconstructed mass distribution. SN image positions 
(and perhaps flux ratios, if microlensing and differential extinction effects are negligible) also 
provide local constraints on the potential, further improving the accuracy of a time-delay Hq 
value. Taken together these considerations suggest that statistical Hq errors on the order of 
a few percent should be attainable if successful observations can be made. And as noted in 
the introduction, the lens model will not be subject to systematic errors associated with a 
global mass sheet degeneracy as long as multiple-image constraints arc available for sources 
at different redshifts. 

To estimate the error that could arise due to time-delay contributions from unmodeled 
galaxy-scale substructure, we note that time delays scale as the square of the characteristic 
deflection of the lens, which in turn scales as the square of the lens velocity dispersion, as in 
(5) and (6). With a time delay of O.lrduster, a cluster velocity dispersion of 1000 km/s, and a 
perturbing galaxy velocity dispersion of 200 km/s, the induced fractional error in Hq would 
be 

SHo ^ Viaxy ^ (200)^ ^ ^ a, (^r^ 

Ho O.lwcr 0.1 X (1000)4 ■ ' ■ ^ > 

Thus substructure is unlikely to be a source of large error in the determination of Hq from 
an observed time delay. An appreciable substructure perturbation to the overall cluster 
potential would also affect lensed SN image positions (which can themselves constrain the 
potential), so it could not go entirely unmodeled by a reconstruction that makes use of all 
lensing constraints. 



5. Predicted Detection Rates 

We wish to quantify detection rates for observable SN image pairings with acceptably 
short time delays. Since type II (core-collapse) SNe will be detected with much higher fre- 
quency than type la's (SuUivan et al. 2000, for example), we focus exclusively on the former. 
Assuming a Salpeter IMF with lower and upper mass cutoffs of 0.1 Mq and 125 Mq respec- 
tively, and that all stars with mass above 8Mq result in core-collapse SNe, Madau (1998) 



derives a conversion factor of 0.0074 between solar masses of star formation and eventual 
number of type II SNe. The comoving star formation rate density at a given redshift can in 
turn be related to an observable quantity such as comoving rest-frame ultraviolet luminosity 
density, which is seen to rise sharply out to the redshifts of ^ ~ 1 that we would hope to 
probe in our search (Madau, Pozetti, & Dickinson 1998). 

Here we consider two versions of the cosmic star formation rate density as a function 
of redshift. First, following Sullivan et al. (2000, "SFH-I"), we use the form of Madau & 
Pozzetti (2000) corrected to a lower mass cutoff of 0.1 Mq and adjusted from an Einstein- 
DeSitter (EdS) universe to our assumed ACDM cosmology in the manner described by Hogg 
(2001). This form incorporates an upward correction to account for extinction (Madau 2000; 
Madau & Pozzetti 2000). We also consider the greatest and least of the three star formation 
rate density functions reported by Lanzetta et al. (2002) based on their determination of the 
star formation rate intensity distribution function from redshifts z — 0-10, again converted 
from EdS to ACDM. These authors do not attempt an extinction correction, and we work 
directly with their "unobscured" star formation rate densities. 

We can combine our assumed SN rate densities with multiple-imaging cross sections 
c(//min, Atmax) computed numerically for the SIE and NIE in an integral over source redshifts 
to obtain an estimate of the detection rates that we might expect: 



dVc/dfldz is the comoving volume per unit solid angle per unit redshift, the angular cross 
section depends upon source redshift and includes a contribution from each independent 
image pairing, psT<s{z) is the rate density of core-collapse SNe per unit comoving volume per 
unit proper time, and the factor of (1 + z)^^ converts from proper time at the source to 
observer time. In an actual calculation this integral must be cut off at some upper limit 
beyond which the supernova rate density becomes utterly unknown. 

To obtain a numerical estimate, we assume a 0"^, = 1000 km/s Icnsing cluster at redshift 
= 0.2 with an isodensity axis ratio q = 0.65 (as mentioned above, with the exagger- 
ated ellipticity we hope to approximate the effects of cluster substructure). We take the 
maximum tolerable observer- frame time delay between multiple images to be 3 hj^ years, 
and investigate a range of differences between limiting detectable apparent magnitude 
and SN II absolute magnitude Msn- We model the SN II spectrum as that of a 10,000-K 
blackbody (Fihppenko 1997), and compute AB- magnitude /c-corrections for observations at 
8140 A (7-band). We consider both an SIE model and an NIE model with a core radius 
of SOh^Q kpc (subtending about 9.1 arcseconds). We normahze the NIE model for zs — 0.7 



SN image pairings 
per cluster 
per obs. year 




(16) 
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as described in Section 3.3, taking this as a typical giant-arc redshift from Williams et al. 
(1999). Figure 3 shows the result of a detection-rate integration out to a redshift of z = 8. 
For visual clarity, we do not plot the NIE rates for the two Lanzetta et al. (2002) star- 
formation histories; they show the same increase by a factor of ~ 2 over the singular case 
as do the Madau & Pozzetti (2000) NIE rates. The results shown are relatively insensitive 
to the choice of upper integration limit, since the star formation rate density of Madau & 
Pozzetti (2000) drops steeply at high redshift, and the A;-correction imposes an effective cut- 
off on any contribution from the enhanced high- 2; star formation reported by Lanzetta et al. 
(2002). 

Although we have considered type 11 SNe as our sources, a rough idea of the detection 
rates for type la SNe image pairings can be obtained from Figure 3 by assuming that type 
la's are ~ 10 times less frequent and ~ 1 magnitude brighter than type IPs. Under these 
assumptions, detection rates for type la events are significantly lower than for type IPs, and 
the /c-correction appropriate to the intrinsically redder colors of type la's will push their 
rates lower still. 



6. Outlook 

The estimated detection rates shown in Figure 3 are far in excess of what one would pre- 
dict based on a circular-cluster lens approximation, and they allow us to gauge the feasibility 
of a cluster monitor program targeted to detect multiply-imaged SNe for use in determining 
Hq. Repeated imaging would be required every month or so, as type II-L (II-P) SNe spend 
~30 (~50) days within one magnitude of maximum light (Doggett & Branch 1985). This 
time will be stretched by (1 + -Zsn), but a precise temporal measurement of peak light would 
be crucial to a time-delay measurement. Taking miim = 24 and Msn = — 18, the most op- 
timistic assumptions predict on the order of 2-3 detectable image pairings per cluster per 
century. A cluster monitor program operating for two to three years would need to image 
~50-100 cluster fields to the required depth in order to have a good chance of detecting a 
few image pairings. Such a program would be feasible with today's 6-8m-class telescopes 
with wide-field cameras given several dark nights per month. The project is better suited to 
a telescope such as the proposed LSST, which would operate in a dedicated survey mode, 
repeatedly imaging large areas of the sky (including strongly lensing clusters) to significant 
depth. Although our predicted image-pair detection rates arc sensitive to assumptions about 
cluster structure and cosmic star formation history, this method of determining Hq certainly 
warrants serious consideration when such observatories as LSST become operational, par- 
ticularly if the current discrepancy between Hq values from lensing and from other methods 
remains unresolved. 
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We stress that in order for this method to be successful, known strong-lensing clusters 
should be targeted specifically. The best constraints on the cluster potential would ideally 
be obtained by follow-up observations with instruments such as ACS of clusters for which 
SN image pairs are actually observed. However, the co-added images from a long-term 
monitor program would be of significant depth, and could also be searched for strongly 
lensed background features. 

The authors wish to thank Paul Schechter for valuable discussion of these topics, Ian 
Small for pointing out the work of Giraud (1992), and the anonymous referee for comments 
that led to a greatly improved manuscript. 
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Fig. 1. — Magnification-limited lensing cross sections. Solid lines are for a g = 0.65 SIE quad. 
The upper, middle, and lower curves correspond to the 2nd/3rd, lst/2nd, and 3rd/4th image 
pairings respectively. Dashed lines are for the corresponding images of a g = 0.65 NIE at the 
"umbilic catastrophe" and with the normalization 6^ = 9e as described in the text. (The 
curve with the reduced maximum cross section is for the 3rd/4th image pairing.) The dotted 
curve is for the two images oi a q — 1 SIS, for comparison. 
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Fig. 3. — Predicted type II SN image pairing detection rates vs. difference between limiting 
observable magnitude and SN absolute magnitude, (/-band, = fOOOlcm/s, q = 0.65, 
= 0.2, maximum observer-frame time delay of Sh^^ years). Solid lines/squares: star 
formation rate density of Madau & Pozzetti (2000). Dashed lines/circles: star formation 
rate density of Lanzetta et al. (2002), upper and lower curves. Open symbols are for an SIE 
model; solid symbols are for an NIE model with a SOh^^ kpc core radius. 



